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Comment on "Liquid Limits: The 
Glass Transition and Liquid-Gas Spinodal 
Boundaries of Metastable Liquids" 

In a recent papeifl, Sastry studies by numerical simula- 
tions the phase diagram of a simple fragile glass-forming 
liquid, the binary mixture with Lennard-Jones (LJ) inter- 
actions introduced by Kob and Andersen, presenting very 
interesting and clear results. Among other things, he ob- 
tains the position of the liquid-glass transition tempera- 
ture as function of density p. He studies both the temper- 
ature To (p) where the diffusion coefficient (of A particles) 
drops to zero according to the Vogel-Fulcher-Tammann- 
Hesse law, and the temperature where the extrapolation 
of data on the configurational entropy Sc{p,T) vanishes, 
Sc being evaluated as the difference between the liq- 
uid entropy and the harmonic amorphous solid entropy. 
These two ways give values remarkably close. 

We apply to the LJ binary mixture, at various density 
values, the analytic approach to-, structural glass ther- 
modynamics recently introducedEI, which allows to ob- 
tain both the liquid-glass transition temperature Tk and 
properties in the glassy phase below Tk doing a first 
principles computation. In one version of this approach, 
which implements analitically the proposal of inherent 
structure dominance at low temperature, one needs to 
compute the liquid free-energy density / and two-point 
correlation functions g{r) at T ~ Tk- 

The theory and different wajya for implementing it are 
discussed in detail elsewhereBij. We only recall that, 
within this framework, at Tk the system undergoes a sec- 
ond order transition corresponding to a one step replica 
simmetry breaking. As one could expect intuitively, it 
is found that in the region T ~ Tk, where the phase 
space can be partitioned in 'valleys', the entropy of the 
liquid corresponds to the entropy of a typical valley Ssoi 
plus the configurational entropy Sc, and the Kauzmann 
temperature is defined by the condition Sc{Tk) = 0. 

We use the harmonic approximation scheme, which 
gives a simple formula for the single valley entropy in 
the liquid phase: 

Ssoi = ^ In (27r e T) - i (Tr In M) , 

where A4 is the 3A^ x 3A^ Hessian (i.e. the matrix of sec- 
ond derivative of the potential with respect to the N d 
particle coordinates) and (•) means the Boltzmann aver- 
age. We are therefore approximating 5*50; by only keeping 
the harmonic vibrational contribution. Further approx- 
imations are needed for computing the spectrum. Here 
we consider fluctuations of the diagonal part of A4 up to 
second orderO, while the whole nondiagonal contribution 
is summed up by means of a kind of chain approximation. 

In the computation of the liquid free energy and of the 
correlation functions in the high density region we use 
the Hypernetted Chain Approximation (HNC). Unfor- 
tunately this approximation overestimates the gas-liquid 
transition temperature and a more refined one is neces- 


sary at low densities. Liquid quantities at p = 1.2 have 
been computedO within a scheme similar to that intro- 
duced by Zerah and Hansen (ZH), which also has the 
advantage of reducing the error in the approximation (of 
order 10% in the HNC case). 


1.4 

Tq from Diffusivity 

i 


(HNC for liquid quantities) 


1.2 

(ZH for iiquid quantities) 

□ ■ 

1 


■ 

0.8 


■ 


■ 

♦ 

0.6 



0.4 




8 ♦ 


0.2 

♦♦♦♦ 


n 




1.1 


1.2 


1.3 
P 


1.4 


1.5 


FIG. 1. The data on To obtained by SastryU compared with 
our analytical results. 

We plot in [Fig. 1] the theoretical predictions on Tk at 
different density values, compared with the data on Tq by 
Sastry. Despite of the error in liquid quantities obtained 
by HNC, there is a quantitative agreement which makes 
evident that the analytic approach works well and it can 
be used as a powerful tool for a better understanding of 
glassy systems. 

To conclude, we stress that, at least in this case, the 
Kauzmann temperature turns out to be directly compa- 
rable with the temperature where the diffusion drops to 
zero, which strongly supports the theoretical expectation 
of an effective ergodicity breaking. 
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